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ABSTRACT 

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to 
investigate the variability of objects in 8.1 deg 2 of the NOAO Deep Wide Field 
Survey Bootes field. We perform a Difference Image Analysis of the four available 
epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 /im bands. 
Out of 474179 analyzed sources, 1.1% meet our standard variability selection 
criteria, that the two light curves are strongly correlated (r > 0.8) and that their 
joint variance (<Ti 2 ) exceeds that for all sources with the same magnitude by 2a. 
We then examine the mid-IR colors of the variable sources and match them with 
X-ray sources from the XBootes survey, radio catalogs, 24 /im selected active 
galactic nucleus (AGN) candidates, and spectroscopically identified AGNs from 
the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, 
most of the variable sources are AGNs (76%), with smaller contributions from 
stars (11%), galaxies (6%), and unclassified objects, although most of the stellar, 
galaxy, and unclassified sources are false positives. For our standard selection 
criteria, 11%— 12% of the mid-IR counterparts to X-ray sources, 24 /im AGN 
candidates, and spectroscopically identified AGNs show variability. The exact 
fractions depend on both the search depth and the selection criteria. For example, 
12% of the 1131 known z > 1 AGNs in the field and 14%-17% of the known AGNs 
with well-measured fluxes in all four IRAC bands meet our standard selection 
criteria. The mid-IR AGN variability can be well described by a single power- 
law structure function with an index of 7 m 0.5 at both 3.6 and 4.5 /im, and an 
amplitude of So — 0.1 mag on rest-frame timescales of 2 years. The variability 
amplitude is higher for shorter rest-frame wavelengths and lower luminosities. 

Subject headings: cosmology: observations — galaxies: active — infrared: galax- 
ies — quasars: general 

1. Introduction 

A substantial fraction of objects in the universe change in brightness with time. This 
apparent variability can be due to physical periodic changes in the objects (e.g., Cepheids, RR 
Lyrae), motion (e.g., eclipsing binaries, microlensing, rotating stars), explosive events (e.g., 
supernovae, novae), and accretion (e.g., cataclysmic variables, active galactic nuclei (AGNs)). 
In a high-latitude extragalactic survey, the most common variable sources are active galactic 
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nucle i (e.g ., Sesar et al. 2007 ). O ptical variability is on e method for identifying quasars^ 



Everil2002l ; iGeha et al. 



20031 ; iRengstorf et al.l I2004J ) , although it is o nly recently that 



(e.g., 

a ful ly quantitative approach to variability selection has been developed ( IKozlowski et al. 
2010h . 



The emission from quasars generally has three components. The UV to near-IR radiation 
is dominated by a hot accretion disk extending from an inner edge of a few gravitational radii 
from the black hole outwar d with, in simple thin disk theory, a temperature profile T oc i?~ 3 / 4 
(IShakura fc Sunyaevlll973l ) . Near the inner edge ; ther e is a corona of hotter gas that produces 
the non-thermal X-ray emission (IHaardt et al.l 1 19941 ). On scales where the temperature is 
below the dust sublimation temperature (~ 2000 K), dust absorbs radiation from the disk 



and reradiates the energy in the mid-IR and far-IR (jBarvainisI Il987l ) . The overall spectrum 
typically shows a minimum near 1 /xm, with the emission from the d isk rising toward th e 
UV and the emission reradiated by dust rising toward the far-IR (e.g., ISanders et al.lll989l ). 
There is an increasing evidence that many physical relations between AGNs, galaxies and 
their large scale cluster ing have to be tak en into account if their formation and evolution is 
to be understood (e.g., iTasse et al.ll2008f ). 



We know a great deal about the optical variability of quasars both from large studies 
of the variability seen in ensembles of sparsely mo nitored quasars and from detailed stud- 



ies o f individual quasars. Ensemble studies (e.g., IVanden Berk et al.l 12004 ; Ide Vries et al. 
20051 ) have shown that variability increases with decreasing optical wavelength, decreasing 
luminosity, and potentially decreasing black-hole mass. The structure function of the en- 
semble variability is a power law with smaller variability amplitudes on short timescales, 
with some evidence for saturation on timesc ales of order a few decades. Until recently, there 
were few studies of individual quasars (e.g., ICutri et al.lll985l ; IClavel et al.lll989l ; iHook et al. 
19941 ). but this has changed dramatically in the last year. The light curves of individual 
quasars are well modeled by a damped random walk, a stochastic process described by the 
amplitu de of the random walk and a damping timescale for returnin g to the mean lumi- 



nosity flKellv et al.l 120091; IKozlowski et al.l I201CH : MacLeod et al.ll2010f ). While preliminary 
indications from Kelly et al.l ( 120091 ). based on ~ 100 quasars, suggest that these two process 



P 
• ( 



parameters are related to the quasar luminosity and black hole mass, iMacLeod et al.l (120101 ) 
used ~ 9000 Sloan Digital Sky Survey (SDSS) Stripe 82 quasars to find a number of clear 
trends. For example, the asymptotic variability on long timescales decreases with increasing 
luminosity and rest-frame wavelength, and is correlated with black hole mass. The timescale 
for returning to the mean luminosity increases with wavelength and also with increasing 



1 We will use words "quasar" and "AGN" interchangeably throughout this paper. 
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black hole mass, but remains constant with redshift and luminosity. 



Far less is known about the near-IR and mid-IR variability of quasars. Since AGN lumi- 
nosities do vary in time at optical wavelengths, we expect to find some variability in the IR, 
but it could be heavily smoothed by averaging the re s ponse over the large scales of either the 
cooler parts of the disk or the absorbing dust. I Glass! ( 12004] ) found that the majority (39/41) 
of the surveyed low-luminosity nearby Seyferts do vary in the near-IR and mid-IR, and that 
the variability is most apparent at longer wavelengths due to diminishing flux from the host 
g alaxy. The most extensive work t o date on near-IR and mid-IR variability of quasars is that 
of iNeugebauer fc Matthews! (Il999l ). These authors studied 25 low-redshift quasars (z ~ 0.1), 
which were observed in five bands (1.3 - 10.6 /jm) for ~ 30 years. The amplitude of variations 
for radio loud (quiet) quasars is ~ 0.3 mag (0.1 mag) on rest-frame time lags (time differ- 
ences between any two epochs) of 2 years at 10.6 /mi. For one quasar, PG 1226+023, they 
measured structure functions in five bands, all showin g ~ 0.1 mag varia t ions f or a rest-frame 
time lag of 2 years. Another example is the work of ISuganuma et al.l ( 120061 ). who found a 
delayed response between the K- and V-band variability in four Seyfert 1 galaxies. These 
authors showed that .fT-band light curves have a fairly tight time lag-luminosity relation, 
r f=s L 1 / 2 , while broad line region (BLR) time lags are more scattered. The interpretation is 
that the i^-band emission is mainly coming from dusty clouds bounding the BLR, and that 
the size scale of these dust clouds is very dependent on the AGN luminosity, since they have 
to be colder than the dust sublimation temperature. 

In this paper, for the first time, we investigate the variability of a truly large num- 
ber of mid-IR sou rces, using the four-epoch Spitzer Deep Wide Field Survey (SDWFS; 
Ashbv et all l2009h o f the Bootes field of the NOAO Deep Wide Field Survey (NDWFS; 
Jannuzi fc Deyill999l ). The two goals are to characterize the mid-IR variability of a large 
sample of AGNs , and to use the variability to identify additional AGNs. While mid-IR 
quasar selection (IStern et al.ll2005l ) works very well when the AGN dominates the luminos- 
ity (e ven for quasars behind dense stellar fields such as LMC/SMC; iKozlowski &: Kochanek 
2009 |), the method begins to fail as the luminosity becomes comparable to that of the host 
(IGorjian et al.ll2008l ; lAssef et al.ll2010bl ). With mid-IR variability, we should be able to iden- 
tify such sources and extend the sample of mid-IR selected quasars to lower luminosities. We 
also make use of th e extensive wavelength coverage of the field, including X-ray (XBootes; 
Murray et al.l 120051 ). the earlier Infrared Array Camer a (IRAC) Shallow Survey that be- 
came the first epo ch of SDWFS (lEisenhardt et al.l 120041) . the 24 fim s urvey of the field (see 
Houck et al.ll2005h . and the FIRST faecker et al.lfl995h and WSRT Jde Vries et all l2002h 
radio surveys. Based on these data, the AGN and Galaxy Evolution Survey (AGES; C. S. 
Kochanek et al. 2010, in preparation) selected galaxies and AGNs at all th ese wavelengths 



for spectroscopic follow up using the Hectospec instrument on the MMT (IFabricant et al. 



-5 - 



20051 ). providing roughly 23,500 redshifts in the field. In Section [2] we describe the data 
analysis, and in Section [3] we describe our variability selection criteria. In Section HI we 
present the results of this study. The mid-IR structure function analysis is presented in 
Section followed by a summary in Section Throughout this paper, we use a standard 
ACDM model with {Q A ,Q M ,Q k ) = (0. 7,0.3,0.0) and h = H n /100 = 0.70. Where needed, we 
use the galaxy and AGN templates of Assef et al. ( 2010bl ) to describe the mid-IR colors of 
various populations and to make any if-corrections or estimates of absolute luminosities. 



2. Data Analysis 



The SDWFS is d escribed in detail by lAshby et al.l (120091 ). The images were taken with 
the IRAC instrument ( Fazio et all2004 ) on the Spitzer Space Telescope, covering wavelengths 
of 3.6, 4.5, 5.8, and 8.0 /im, hereafter referred to as the [3.6], [4.5], [5.8], and [8.0] bands, 
respectively. A single IRAC frame is 256 x 256 pixels. We work with full- field SDWFS mosaics 
constructed from the ~20,000 IRAC images/epoch. Each approximately 13600 x 16700 pixel 
mosaic covers 10.2 deg 2 , centered at (R.A., decl.) = (14:31: 27, +34:10:43), with a p ixel scale 
of 0.84 arcsec. Counting the original IRAC Shallow Survey (jEisenhardt et al.ll2004j ) as epoch 
1, the SDWFS data span a period of four years between 2004 and 2008, where the intervals 
between the epochs are approximately 3.5 years, 0.5 years, and 1 month (observations were 
done on 2004 January 10-14, 2007 August 8-13, 2008 February 2-6, and 2008 March 6-10). 
For each individual observation of a patch of the sky (256 x 256 pixels), there is only a ~ 40 
s offset between the channel 1/3 and 2/4 observations. For a given SDWFS epoch, each 
patch of sky was observed ~ 3 times over a 1-2 day period. Each SDWFS epoch of the 
entire 10 deg 2 area took ~ 5 days, and we treat it as a single image with a single Julian 
Date. This does not affect our final conclusions, as we show later that there is little mid-IR 
AGN variability on these timescales. 

To study the variab ility of objects in the SDWFS field, w e used the Difference Image 
Analysis (DIA) method (jAlard fc Luptonlll998l : Iwozniaklboool ). While the variability itself 
is the subject of this paper, the DIA method also proved to be a useful diagnostic tool for 
constructing artifact-free mosaics (see lAshby et al.l 120091 ) . The basic idea behind the DIA 
method is to match one image, called the template or reference image, both astrometrically 
and photometrically to another image of the same area of the sky. Once they are astromet- 
rically aligned, the template image (which is chosen to be of better quality) is convolved 
with a kernel function to match the point spread function (PSF) of the analyzed image. The 
two images are then subtracted from each other (Figure [I]) to leave only variable objects, 
noise, and systematic problems (cosmic rays, artifacts, etc.). After some experimentation, 
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we concluded that the [5.8] and [8.0] maps were too shallow to perform DIA and obtain 
reliable light curves, so our variability analysis is restricted to the [3.6] and [4.5] bands. 

Spitzer PSFs have a "triangular" shape with extended diffraction-like patterns (more 
precisely they are "multiplexer bleed" or "muxbleed," "column pull-down," and "banding''^), 
and the PSF rotates between the epochs. DIA works best on simple quasi-Gaussian PSFs, 
so we pre- filtered the images to have a simpler PSF structure (Figured]). For each region, we 
scaled and median combined the brightest stars to build an empirical model of the PSF, B. 
We then linearly filtered the image / to create a filtered image defined by the combination 
of Fourier transforms If = GI/B, that transforms the PSF from B to G, where G is a 
Gaussian broad enough to avoid noise amplification (i.e., G — > faster than 5—^0 for small 
scales/large wave numbers). The filtering process conserves fluxes and was carried out on 
a 3 x 4 grid of 4040 x 4040 pixel regions that overlap by 20 pixels along their edges. To 
save some computation time, each of these 4040 x 4040 regions was then cut into smaller 
2040 x 2040 pixel sub-images (again each with 20 pixels of margin) for the DIA analysis. 
We combined all four epochs to build the te mplate images and these templates were scaled 



and subtracted from each epoch using DIA (jWozniak! 120001 ) . The brightest stars in the field 



produce artifacts due to subtraction residuals, particularly from the "diffraction" spikes in the 
PSF. We therefor e masked a region around all stars with K < 11 mag in the 2MASS catalog 



of point sources ( ICutri et al.l 120031 ). The masking radius was 260(1 — K/ll) arcsec, which 
is large enough to include the spikes. We used the Two Micron All Sky Survey (2MASS) 
catalogs because these bright stars were saturated and not included in the SDWFS catalogs. 
The masked regions constitute ~ 3% of the area. 

The [3.6], [4.5], [5.8], and [8. 0] Vega magnitude s, and the coordinates (J2000.0) are 



taken from the SDWFS catalogs (lAshby et al.l 120091 ). Throughout this paper, we use 4 



arcsec diameter aperture magnitudes with aperture corrections to 24 arcsec. The original 
3.6 /zm full SDWFS catalog contains 677,522 objects (catalog DR1.1), but we cropped it to 
515,055 objects (hereafter called the "cropped sample") as variations in coverage between 
epochs reduced the usable area for variability analysis by about 20%, to 8.1 deg 2 . Using 
DIA's PSF fitting routines, we could construct light curves (at least three epochs) for 492,420 
objects simultaneously in both [3.6] and [4.5]. The final number of objects used in this 
analysis after masking regions around bright stars is 474,179. While we include the results 
for all objects from the cropped sample in Table [lj in the analysis that follows we consider 
only the sources with at least three epochs. Table [2] provides the DIA light curves for all 
SDWFS sources from the cropped sample. 



2 For details see the IRAC Data Handbook at 



http : //ssc . spitzer . caltech . edu/irac/dh/ 
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In order to understand the nature of the variable sources, we matched the SDWFS 
sources to several of the available data sets for the Bootes field, focusing on information 
related to AGN activity. Thus, we matched to the X-ray and radio catalogs, to the subset 
of 24 /iin sources with high probabilities of being AGNs, and to the AGES spectroscopic 
data. We used the X-ray catalogs from the XBootes Chandra survey covering 9.3 deg 2 of the 



Bootes field in the .5-7 keV band (IMurrav et al.1120051: iKenter et al 120051 : iBrand et al.ll2006 
Gorjian et al. 2008 ). We matched a version of the Murray et al. ( 20051 ) catalogs extended 



downward in flux from 3293 sources with 4 or more counts to 4642 sources with 2 or more 
counts. While using 2 count (about 4x 10~ 15 erg cm -2 s™ 1 ) X-ray sources sounds peculiar, the 
backgrounds in XBootes are so low that almost all such sources with mid-IR identifications 
will be real, particularly if the match is to one of the relati vely rare mid-IR variable sources. 
Using the Bayesian matching method of IBrand et al.l (120061 ) . we identified 3363 X-ray sources 
with at least t wo epochs of SDW FS data. We used th e same procedures t o match to the 
20 cm FIRST faecker et al.lll995h and 1.4 GHz WSRT (jde Vries et al.ll2002h radio catalogs. 
Restricting the matching to likely radio point sources in the FIRST and WSRT surveys, we 
found 322 and 2118 matches, respectively, to SDWFS sources with at least two epochs of 
data. We know from AGES that most 24 urn, optical point sources satisfying the criterion 



that I > 18 — 2.5 log(-F24/mJy) are quasars (iBrown et al.ll2006l : lEckart et al.ll2010l ; lAssef et al. 



2010aJ), where the optical criterion is used to eliminate stars. Finally, we matched to the 



AGES redshift survey, finding 19,885 (18,288, 1387, 497) SDWFS sources (wit h redshifts 
z > , 1,2). The modified SDSS pipeline used for the AGES redshifts (HSRED; ICool et al. 



20081 ) provides an estimate of the best fit spectral template used for the fits that can be used 



to classify objects. More generally, we have all the galaxy and AGN target selection codes 
used by AGES to select and prioritize spectroscopic targets. 



3. Variability 

For each light curve i (from Table [2]) with at least three epochs in both the [3.6] and 
[4.5] bands, we computed the standard deviation t>([X]) J for that light curve as a function of 
(Vega) apparent magnitude m([X])*, as shown in Figure [2j We also calculated the variability 
covariance C{ 2 and Pearson's correlation r l coefficients between the two channels, where 

N \ l / 2 

1 " epochs 



v([X]T = - 1 — jr (m([X])}-(m([X]y>) 2 , (1) 

- -^epochs 

C\ 2 = — - (™([3.6]y» (m([4.5])} - (m([4.5]y» , (2) 



* epochs j~l 
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ni 

v{[3.6]yv{[4.5}y , y J 

and (m([X])*) is the average magnitude for the ith light curve in channel X. The correlation 
coefficient is limited to the range — 1 < r < 1, where r = 1 (—1) means that the [3.6] and 
[4.5] light curves are perfectly correlated (anti-correlated). The distribution of t>([X])* is 
dominated by contributions from noise (and any systematic errors) rather than variability, 
so we need to define a threshold for selecting candidate variables from these distributions. 
Let f m ([X]) be the median of the t>([X]) 1 as a function of magnitude and cr([X]) be the 
standard deviation of the v {[X]) 1 values around this median. Sources are more likely to be 
true variables as their distance, (v ([X]) 1 — v m ([X])) /a([X]) > 0, from the median increases, 
and if the [3.6] and [4.5] variabilities are strongly correlated. A reasonable metric for the 
joint significance of the variances in the light curves is 

^([4.5])--^([4.5]) \ 2 \ 1/2 

-([4.5]) ) ) ' W 

which quantifies the degree to which the object deviates from the median variances in both 
bands. Objects with positive excess variances (large a 12 and v > v m ) and strong correlations 
r ~ 1 between the two bands are the best candidates for true variables. 

Figure [3] shows the distribution of correlation strengths r for sources with 012 exceeding 
a range of thresholds. We see that there is a peak for highly correlated sources r ~ 1 whose 
strength increases as the variability significance 012 increases. This demonstrates that the 
variability in the two bands is tightly correlated for real variables and that we are likely to 
produce a relatively clean sample of variables (the false positive rate is Q%-7% for AGNs 
in the Photometry Group 3 defined below) if we restrict ourselves to r > 0.8 and a 12 > 2. 
Variability selection is always a tradeoff between completeness and contamination, so we 
will consider sources in three different variability levels, all restricted to have correlations 
r > 0.8 between the bands. To avoid confusion, the Variability Levels correspond directly to 
the values of Oyi- Variability Levels 2, 3 and 4 correspond to sources with excess variances 
and significances exceeding 12 > 2, 3 and 4, respectively, with Variability Level 2 being the 
most complete and the most contaminated, and Variability Level 4 being the least complete 
and least contaminated. 

We estimate the fraction of false positives by assuming that the density of uncorrelated 
sources (r < 0.5) in Figure [3] represents the density of false positives amidst the highly 
correlated (r > 0.8) sources. If the number of identified variables is iV r> o.8, that is sources 
with > 2 (or 3, 4) and r > 0.8, and the number of sources with 012 > 2 (or 3, 4) and 
r < 0.5 is N r<0 _5, then we estimate that the number of false positives is F = (0.2/1.5)A^ <0 5 , 



j / ^([3.6])--^ m ([3.6]) \ 2 | / 
I V <x([3.6]) J V 
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where 0.2/1.5 is the ratio of the intervals r > 0.8 and — 1 < r < 0.5. The fraction of the 
identified variables that are false positives is then F/N r>0 . 8 . 

Our ability to characterize sources also depends on the magnitudes of the object, because 
this determines how well we can measure the mid-IR colors. We divide sources into three 
photometry groups based on their mid-IR magnitudes. 



• Photometry Group 1 consists of all objects that after masking have three or four epochs 
of data, with no constraints on either their magnitudes or magnitude uncertainties. 
There are 474,179 Photometry Group 1 objects after masking, of which 5107, 2208, 
and 1071 meet the Variability Level 2, 3, and 4 criteria. 

• Photometry Group 2 is the subset of the Photometry Group 1 sources for which we 
can measure the [3.6] — [4.5] color. These sources must satisfy [3.6] < 19.7 mag and 
[4.5] < 19.3 mag with uncertainties in these bands smaller than 0.1 mag. There are 
213,594 Photometry Group 2 objects, of which 1557, 668, and 359 meet the Variability 
Level 2, 3, and 4 criteria. 

• Photometry Group 3 is the subset of Photometry Group 2 objects for which we can 
also measure the [5.8] — [8.0] color. These sources are required to have uncertainties 
smaller than < 0.2 mag for the 5.8 and 8.0 /j,m bands. There are 39,522 Photometry 
Group 3 objects, of which 775, 428, and 260 meet the Variability Level 2, 3, and 4 
criteria. 



Table [TJ provides these significance statistics for the SDWFS sources. Note, however, that 
our catalogs are limited to the objects detected on the template image - any transient which 
does not produce a detectable source in the template image, which is just the average of four 
epochs, will not be identified. The most variable, highly correlated (r > 0.8) sources in the 
field are AGNs, peaking at cr 12 w 20. 

We will make extensive use of the [3.6] — [4.5] versus [5.8] — [8.0] mid-IR color-color dis- 
tribution to characterize the sources. Figure B] shows th e distr ibution of all Photometry 
Group 3 SDWFS sources in this space. lEisenhardt et al.l ( 120041 ) suggested that for bright 
sources the r ed plume in 3.6] — [4.5] color was likely to be dominated by AGNs. This was 
confirmed bv lStern et al.l (20051) b ased on the first year of spectroscopic observations by the 



AGES survey. IStern et al.l (l2005f ) defined a selection region, the "AGN wedge" , for mid-IR 
AGN selection that is also shown in Figure HI For these Vega colors, stars lie at color 
unless they have dusty winds o r disks to produce a mid-IR excess. Using the spectral energy 
distribution (SED) models of lAssef et al.l (J2010b|), which are derived by fitting the SEDs 
of Bootes sources, we can illustrate the locations of various extragalactic populations as a 
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function of redshift. Early-type galaxies at low redshifts have colors similar to stars and 
then become redder in [3.6] — [4.5] for z > 1 at roughly fixed [5.8] — [8.0], tracking along the 
left edge of the AGN wedge. They become red in [3.6] — [4.5] due to the flux "peak" near 
rest-frame wavelength of 1.6 /im, which passes through the 4.5 /zm band at a redshift of 1.5- 
2.1. Late-type galaxies at low redshift have red [5.8] — [8.0] colors at low redshift due to the 
polycyclic aromatic hydrocarbon (PAH) emission in the 8 jum band, then move horizontally 
to bluer [5.8] — [8.0] along the bottom of the AGN wedge as the feature is redshifted out of the 
8 /xm band, and finally follow the early-type galaxies up along the left edge of the selection 
region. Pure AGNs remain inside the AGN wedge except near z ~ 4.5 where they make an 
excursion out of the AGN wedge when the Ha line lies in the 3.6 /im band. If we use an 
equal combination of an AGN and a host galaxy, defined by the total emission from 0.1 to 
30 fim, the combined SED is somewhat bluer. With further increases in the host luminosity 
relative to the AGN, the source drops below the bottom of the AGN wedge. 

In fac t, the distribution of the IRAC co lors for the XBootes X-ray sourc es does extend 



below the IStern et al.l ( 120051 ) AGN wedge ( jGorjian et al.l 120081 ). Based on iGorjian et al. 



( 120081 ) . we define a "modified AGN wedge" that encompasses most X-ray sources. This 
modified AGN wedge is defined by [3.6] - [4.5] < 2.5 x ([5.8] - [8.0]) - 0.6 mag (left side), 
[3.6] - [4.5] > 2.5 x ([5.8] - [8.0]) - 3.5 mag (right side), and [3.6] - [4.5] > 0.2 x ([5.8] - 
[8.0]) + 0.08 mag (bottom). Because we are identifying AGN candidates based on their mid- 
IR yariability, we sho uld be able to identify them in this [3.6] — [4.5] blueward extension of 
the lStern et al.l (120051 ) selection region despite any host contamination. The caveat, however, 
is that this blueward extension is dominated by false positives (galaxies) if used to select 
AGNs based solely on their mid-IR colors without additional information such as variability 
or X-ray emission. We also define a "stellar" box which should contain most normal stars. 
The edges of this box are defined by —0.15 mag < [3.6] — [4.5] < 0.15 mag and —0.15 mag 
< [5.8] — [8.0] < 0.15 mag. The coldest stars, reaching into the brown dwarf regime, become 
red in [3.6] — [4.5] due to methane absorption in the 3.6 /x m band, but have [5.8] — [8.0] colors 



near zero and thus stay to the left of the AGN wedge (IStern et al.l 120071 ; lEisenhardt et al. 
2010h . 



Beyond these mid-IR color diagnostics, we also use the extensive multi-wavelength and 
spectroscopic data for the field to characterize the sources. We know that the variable 
sources must be dominated by AGNs and a small fraction of stars, so the basic check on our 
selection criteria is that our selection limits on a±2 and r preferentially select sources with 
AGN signatures, without introducing large numbers of contaminating sources. 
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4. Results 

We will now characterize the variable sources using three approaches. First, we simply 
examine the distribution of objects in mid-IR color. Second, we examine the variability of 
AGN candidates identified at other wavelengths. Finally, we consider the spectral properties 
of the objects with known redshifts. 

Figure |5] and Table [3] characterize the distribution of objects in mid-IR color. The 
Photometry Group 3, Variability Level 2 candidates are concentrate d in either the modi fied 



sec 



AGN wedge, with a distribution very similar to that of X-ray sources (jGorjian et al.ll2008 
Figure H]), or in the stellar box. While only 0.5% of the Photometry Group 3 objects outside 
these two regions are Variability Level 2 sources, 7.4%, 4.6%, and 1.1% of the sources in the 
old AGN, modified AGN, and stellar regions are Variability Level 2 sources. The estimated 
false positive rates are 44%, 6%, 7%, and 58% for these four regions (outside, old AGN wedge, 
modified AGN wedge, stellar box) in the color-color space, respectively. This simply means 
that, for example, of the 7.4% variable sources in the old AGN wedge 6% are false positives, 
so in fact 7.4%x (1-0.06) = 7.0% of all objects are real variables and 7.4%x0.06 = 0.4% of all 
objects are falsely identified as variable. Increasing the variability significance level reduces 
the variable fractions roughly in the same proportion, suggesting that even the Variability 
Level 2 candidates are dominated by true variables. The distribution also validates our 
assumption that the variability is dominated by AGN - 75% of the candidates lie in the 
modified AGN region for Variability Level 4. The selection remains pure even for Variability 
Level 2 because the requirement that the [3.6] and [4.5] light curves are tightly correlated 
(r > 0.8) eliminates most of the contamination even when we only require Oyi > 2 (see 
Figure [3]). Most variable sources in the modified AGN wedge are known to be AGNs based 
on other evidence. Only 16%, 10%, and 7% of the Photometry Group 3, Variability Level 
2, 3, and 4 sources lack photometric or spectroscopic signs of AGN activity. Most of these, 
however, are likely to be false positives given our estimated false positive rates of 7%, 4%, 
and 3%. For the Photometry Group 2 objects, where we lack accurate [5.8] — [8.0] color, we 
can only investigate the distribution in [3.6] — [4.5] color, and the side panel to Figure shows 
the fraction of variable sources as a function of this color. Not surprisingly, the fraction is 
high (~ 10%) when the color is red, corresponding to AGNs, and then drops dramatically for 
colors consistent with normal galaxies. The fraction of variables in the other two photometry 
groups are dramatically lower than for Photometry Group 3 for two reasons. First, a source 
must be intrinsically more variable in order to meet any variability criterion because the 
intrinsic measurement errors are larger. Second, with less well-measured colors, objects are 
scattered from their true location in the color-color space, and the differential scattering of 
galaxies into the AGN regions and AGNs out of these regions dilutes the selection criterion. 
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The 86 (Photometry Group 3, Variability Level 2) variable objects inside the stellar 
box are red, with optical colors of R — I > 1 (i.e., V — I > 2), so they must be dominated 
by red giants or main-sequence stars. The red and asy mptotic giant branch stars include 
the OGLE Small Amplitude Red Giants (OSARGs; e.g.. ISoszvhski et al.l 120041 ) with /-band 
amplitudes of ~ 0.1 mag and perio ds between 10 and 1000 days, the Long Period Variables 
(LPVs; e.g., ISoszvhski et al.l 120091 ) with amplitudes of up to ~ 3 mag and pe riods 10-1000 
days, and the Long Secondary Period variables (LSPs; e.g., ISoszyhskil 120071 ) with periods 
30-1000 days. However, given their absolute magnitudes of Mj ~ —4 mag, they would have 
to be at a distance of ~ 200 kpc to match the observed R, I, and IRAC magnitudes. Since 
this scenario seems unlikely, the variable objects insi de the stellar box could be the main- 



sequence flaring M dwarfs (e.g., iKowalski et al.ll2009f) with absolut e magnitudes of Mi « 7 
mag, and distances of up to 1 kpc. Based on lKowalski et al.l ( 120091 ). we expect to find ~ 10 



M dwarf flares in our data. The R — I and I — [3.6] colors of the variable object s inside the 
stellar box are inconsistent with the galaxy templates from lAssef et al.l (J2010bJ) . We have 
not pursued these sources in any detail because they are almost certainly quotidian variable 
stars and false positives. The estimated false positive rate of 58% is very high. 

Next, we examine the variability of sources whose X-ray, MIPS 24 /im or radio fluxes 
indicate that the source is an AGN. The color distributions of these sources are shown in 
Figure [6] and the variability statistics are summarized in Table HJ Here we focus on Photom- 
etry Group 1 since we are not dependent on having mid-IR colors to interpret the results. 
First, we see that a significant fraction of the X-ray sources show mid-IR variability, with 12% 
(= 12.1% x (1 — 0.04) accounting for false positives) being variable even with no limits on the 
mid-IR photometric uncertainties. The higher fractions in Photometry Groups 2 and 3 are 
largely due to our improved ability to detect a given level of variability as we require smaller 
photometric uncertainties. The MIPS quasar candidates, which are found to be either AGNs 
(~ 80%) or optically unresolved star forming galaxies upon spectroscopic examination, also 
show high variability rates, with 11% (= 11.7% x (1 — 0.05)) of the Photometry Group 1 
objects showing variability at Variability Level 2. The variable fraction again increases as 
we go to the brighter photometric groups. The two sets of radio sources show si gnific antly 
lower levels of variabil ity. This is consistent with the results of iHickox et al.l ( 120091 ) and 
Griffith fc Sterol (120101 ) who show that, as compared to X-ray and mid-IR selected AGNs, 
radio-selected AGNs have lower Eddington ratios and are much less likely to be AGN domi- 
nated. Thus, variability in radio-selected AGNs will be more challenging to discern relative 
to the photometric noise due to the host galaxy. 

Finally, we compare the variable sources with the template types used to estimate the 
source redshifts in the AGES survey, as shown in Figure [7] and summarized in Table |5j 
This is somewhat redundant with the previous two tests because the objects targeted as 
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AGNs are chosen to be in the IRAC AGN wedge, to show X-ray or radio emission, or to be 
optical point sources with non-stellar 24 /zm fluxes. The success rate for the spectroscopic 
identification of candidate AGNs in AGES is quite high. The sample of stars with spectra 
is very heterogeneous. Roughly half of the stars are F and G stars observed by AGES as 
spectrophotometric calibrators. Most of the remainder are K and M stars which were flagged 
as candidate AGNs either because of their X-ray fluxes or the presence of a mid-IR excess. 
We again see that the fraction of variable sources is highest for the AGNs, with roughly 11% 
meeting the Level 2 variability criterion. 

In summary, small relatively shallow mid-IR variability surveys like SDWFS are a work- 
able but not a highly efficient means of identifying AGNs. Roughly speaking, the efficiencies 
are of order 10% for variability selection criterion providing relatively low false positive rates. 
It d oes, however , overc ome the weaknesses of mid-IR color selection to identify AGNs outside 



the IStern et al.l (120051 ) selection region. With difference imaging, the host galaxy flux that 
moves the source outside this color region ceases to be a problem other than through its 
effects on the photometric noise. 



5. Mid-IR Structure Functions 

The structure function is a measure of the expected changes in an object's flux as a 
function of the time between the measurements. The time lag r is defined as the time 
difference between any two data points in a light curve, where the observed time lag is 
t = \ti — tj\ and the rest-frame time lag is r = |tj — tj\/(l + z). For each time lag, one 
wants to find all pairs of points in a light curve and calculate the variance of the magnitude 
differences, Am = — rrij. Such a variance can also be calculated by averaging over many 
objects of the same class (i.e., AGNs) for each time lag between any two epochs, as we do 
in this paper. 

Given the depth of the AGE S redshift survey, all z > 1 objects are spectroscopically 



confirmed AGNs (see Figure 2 in IStern et al.l 120051 ). We measured the mid-IR structure 



function of z > 1 Photometry Group 3 AGNs to fixed magnitude limits of [3.6] < 16, 17, 
and 18 mag. Only ~ 10% of these objects were actually detected as significantly variable 
sources (see Table [5]), so most of the observed variance in the magnitudes is at the noise 
level. This is a rather different regime from most optical studies, where the variance is 
dominated by variability, and makes it critically important to correctly characterize the 
contribution to the variance from the measurement errors. We do this by matching each 
AGN with four unresolved galaxies of nearly identical observed flux. We chose these galaxies 
to be Photometry Group 2 and 3 objects with colors [3.6] — [4.5] < 0.1 mag and 0.5 mag 
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< [5.8] — [8.0] < 1.5 mag, thus avoiding AGNs. We then compute the structure function as 



T 



1/2 



(5) 



In addition to the simple calculation over the 160, 690, and 1100 AGNs with [3.6] < 16, 
17, and 18 mag, we also calculated the structure functions for 1000 bootstrap resamplings 
of both the AGN and galaxy lists. These give consistent results, and we generally report the 
median result of the bootstrap samples and the region encompassing 68.3% of the samples 
about the median. The small values at short lags are very sensitive to any systematic 
problems in the noise estimate, and while this empirical approach to correctly characterizing 
the noise appears to work well, our estimates of the amplitude of the structure function 
on the longer timescales are probably more robust than our estimate of the slope of the 
structure function. 

Figure [8] shows the resulting observed and rest-frame structure functions for the [3.6] 
and [4.5] bands as a function of magnitude limit. We generally see more variability on longer 
timescales. We fit our structure functions with a power-law 

s(r)=s (^y (6) 

using r = 4 years (2 years) for the observed (rest-frame) estimates. Table |6] presents the 
resulting fits. The [3.6] and [4.5] structure functions are consistent with each other, but also 
significantly steeper than the z-band structure function from an ensemble of SDSS quasars 



(IVanden Berk et al.ll2004f ). with 7 = 0.303±0.035, and higher amplitude, with S = 0.20 mag 



for r = 2 years. 

The steeper mid-IR slope could be explained by the emission being dominated by larger 
physical scales than the optical emission. At these redshifts, the observed frame mid-IR 
corresponds to the rest-frame near-IR, where the emission has contributions from both the 
disk and hot dust. If we focus on the more compact disk emission, the characteristic scale 
of an accretion disk for emission at wavelength A, defined by the point where the disk 
temperature matches the photon energy, is 

jl6 f Kcst\^ 3 ( Mbh \ 2 ^ / L X 



R 1Q i* ^ cm, (7) 

where L/Le is the Eddington ratio and L = rjMc 2 defines the radiative efficiency 77. Taking 
L/Le = 1/3, rj = 0.1, and z — 1, the optical emission for Mbh = 10 9 M Q comes from a 
region only a light day (0.003 light years) in size, while the mid-IR scale is approximately 
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10 light days (0.03 light years). Thus, our shortest lags correspond to the fastest possible 
response time given the size of t he emission region since the typical AGES quasar has Mbh — 
1O 8 5 M fcollmeier et alil2006h . The orbital, thermal a nd viscous timescales are even longer, 
and Kelly et al. (2009, also see iKozlowski et al.ll2010l ) find that the optical variability can 
be characterized by an amplitude and a timescale that might be associated with thermal 
timescales. In any case, the mid-IR emission is coming from scales so large that it is hard to 
respond to any driving process on the timescale of our shortest observed lags. This should 
suppress variability on these short timescales compared to the optical emission, which can 
respond more easily, leading to a steeper mid-IR structure function. 

Next we subdivided the sample into two redshift bins, 1 < z < 2 and 2 < z < 3, and 
then divided the objects in these bins at the median value of their absolute magnitudes. 
Figure M shows the structure functions of these subdivided samples as compared to the best 
fits to the full samples. At least for the well-determined values at longer lags, there is a 
pattern that shorter rest-frame wavelengths vary more than longer rest-frame wavelengths. 
Also for 1 < z < 2 AGNs, the lower luminosity systems vary more than the higher luminosity 
systems. This is consistent with earlier findings in the optical part of the sp ectrum (e.g., 
Hook et al.1 1 1994 lHawkinsI I2OO2I : IVanden Berk et all booi Ide Vries et all l2005h . In general, 
the trend of more variability with longer lags also holds, although there are some exceptions 
in the [4.5] results. 

We also subdivided the data into eight smaller redshift-absolute magnitude bins and 
perf ormed a global fit that includes a dependence on wavelength in the following form (see 
also llvezic et al.l 120041 ) 



-D 



S(t) =A[1 + B((Mj)+ 25)] (-^—) ° (^-) 

\ 2 years / \ 1 fim J 

where A, B, C, and D are parameters of the fit, (Mj) is the ave rage absolute J-ban d mag- 
nitude in a bin (the rest-frame Mj magnitudes are taken from lAssef et al.l (j2010al )). r r f is 
the rest-frame time lag, and (A) is the average rest-frame wavelength in each bin. We either 
fit (1) all the data or (2) exclude the short time-lags (r < 0.1 year) as the structure func- 
tion contains little or no signal in this regime. For case (1) we found A(l) = 0.13 ± 0.01, 
B(l) = 0.19 ± 0.03, C(l) = 0.41 ± 0.05, D(l) = -0.04 ± 0.12, while for case (2) we found 
A(2) = 0.13 ± 0.01, B(2) = 0.19 ± 0.03, C(2) = 0.46 ± 0.06, D(2) = -0.02 ± 0.13. We find 
no detectable dependence of the structure function on wavelength, even though the data 
cover the redshift range of z — 1-3, the absolute magnitude range of Mj = —24.5 to —26.5 
mag, and the rest-fr ame wavelengths A r f = 1-2 //m. As compared to optical studies (e.g., 
MacLeod et al.ll2010l ). the mid-IR structure function appears to have stronger dependence on 
the luminosity (parameter B), and no dependence on wavelength (parameter D). Note that 
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the large uncertainty in D (ap = 0.13) can potentially hide, within the 2-3<r, the wavelength 
dependence of D ~ 0.3 observed in optical studies. Another possibility is that the mid-IR 
wavelength dependence is much weaker than that observed in optical. 

Table [7] shows the variability amplitudes at the observed time lag of 4 years for various 
classes of objects from Figure 0] assuming a fixed slope of r 1 / 2 . Objects with z > 1 are 
AGNs, and show uniform and consistent variability amplitudes of ~ 0.09 mag as a function 
of magnitude in both channels for fixed 7 = 0.5. The X-ray selected sources inside the 
modified AGN wedge are expected to be AGNs, and not surprisingly, they also show similar 
variability amplitudes of ~ 0.10 mag to the z > 1 objects. The mid-IR counterparts of 
the 24 //m (MIPS) selected sources inside the modified AGN wedge, again, have almost 
identical variability amplitudes. The radio selected objects, however, have a drastically 
different distribution in their mid-IR colors (top right panel in Figure H|) from that of AGNs, 
suggesting that a large fraction of them may not be AGNs. The structure functions based 
on the radio-selected objects inside the modified AGN wedge may be driven by non-AGN 
objects, and we find somewhat smaller variability amplitudes of So = 0.05 ± 0.01 mag and 
So = 0.06 ± 0.02 mag in [3.6] and [4.5], respectively. 



6. Increasing Completeness with Future Surveys 

Given the structure function, we can now estimate the expected statistical properties 
of the AGNs in the current survey or any extension by adding new epochs. The structure 
function S(r) is related to the correlation function by C{r) = — S(t), where is the 
variance as r — > 00. This makes it straightforward to compute the mean x 2 of an AGN 
(x 2 ) = ((rrii — (m)) 2 /cr 2 since (miUij) = a 2 5ij + C(|tj — tj\), where a is the measurement 
uncertainty. Because we fit for the mean, (m), the exact value of Voo is irrelevant. For our 
observed frame power-law structure function, the result for the SDWFS data and cadence is 
(x 2 ) = 3 + 1.6S 2 /cr 2 , where S = S (t / 'r )°' 52 , S ~ 0.1 mag, and a is the typical photometric 
uncertainty. The constant 3 is simply the number of degrees of freedom after fitting the 
mean, so a non- variable source would have (x 2 ) = 3. This expression makes it clear why our 
AGN completeness is so low and drops rapidly with magnitude (increasing a, see Figure [2]). 

We can, in fact, simulate the completeness. We can decompose the correlation matrix 
Cij = C(\ti — tj\) into eigenvectors Vj and eigenvalues A 2 . If G(x) is a Gaussian deviate of 
dispersion x, then random realizations of light curves consistent with the structure function 
are 

m = J2^G(k) + G(a) (9) 
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(see iRybicki fc Presall992l ). If we assume the variability in both bands is identical but the 
noise is independent, then we can compute the fraction of such AGNs we would select as 
a function of the noise a using Monte Carlo simulations. Figure [TU] shows the approximate 
completeness as a function of magnitude for sources with r > 0.8 and either cr 12 > 2 or 4. 
We use the observed structure function and the normalized measurement uncertainties to 
track the median values for [3.6] in Figure [2j For simplicity, the [4.5] band is treated as 
being identical to the [3.6] band. As might be expected from the structure of Figure EJ the 
completeness is high when the measurement errors are small compared to the amplitude of 
the structure function and then drops rapidly for fainter sources. 

In general, the statistics of these simulated light curves agree with observations, indicat- 
ing that our estimates of the structure function are correct. However, they also show that we 
are "pushing the envelope" of variability, because for most AGNs the variability amplitudes 
are merging into the measurement uncertainties, as can be seen in Figure [2j 



7. Summary 

In this paper, we performed DIA photometry on the four epochs of SDWFS data. The 
common area of the Spitzer mosaics covers 8.1 deg 2 of the NDWFS Bootes field, and contains 
474, 179 mostly extragalactic objects. We concentrate on the deeper 3.6 and 4.5 jum channels, 
and provide variability catalogs and light curves. These variability data are cross-correlated 
with the AGES redshift survey, X-ray and radio data, and 24 /iin AGN candidates from the 
same area of the sky to study the variability of various classes of objects. 

We study the mid-IR variability of objects based on the significance level of their joint 
variability (cr 12 ) in both the [3.6] and [4.5] channels focusing on the sources that also show 
a high correlation (r > 0.8) between the two channels. Our standard selection criteria 
with <j 12 > 2 identify ~ 5100 objects as variables, which constitute 1.1% of all detected 
sources in the field. We find that the majority (76%) of the mid-IR variable sources in the 
NDWFS Bootes field are AGNs, where we define AGNs to be objects that satisfy any of 
the following criteria: (1) lie inside the modified AGN wedge, (2) are X-ray, 24 /im or radio 
AGN candidates, (3) have z > 1, or (4) have an AGN spectroscopic template in any region. 

Amongst all the extragalactic classes of objects, AGNs are the primary sources able to 
produce significant changes in their mid-IR luminosity over a several year timescale. Thus, 
mid-IR variability may be used as a new tool for selecting AGNs, especially if combined 
with other methods. However, in our relatively shallow few-epoch survey, only ~ 15% of 
AGNs are sufficiently variable to be detected; therefore, variability-selected samples of AGNs 
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are highly incomplete. This incompleteness is caused by the typical variability amplitude 
being comparable to the survey photometric errors, and is not due to the small number of 
available epochs. For our standard criteria, 14% of X-ray sources, 17% of 24 /zm selected AGN 
candidates, and 15% of spectroscopically confirmed AGNs are detected by their variability. 



The variable AGNs primarily occupy the mid-IR AGN selection wedge of IStern et al. 



( 120051 ). with an extension to bluer [3.6] — [4 .5] colors where the host galaxy domin ates the 
mid-IR colors of low - luminosity AGNs (see iGorjian et al.l l2008t lAssef et al.l l2010b[ ). If we 
redefine the original IStern et al.l ( 120051 ) selection wedge to encompass X-ray or variability- 
selected AGNs, we find that variable objects bluer than [3.6] — [4.5] > 0.3 mag are likely to be 
AGNs. This opens a new window for the Spitzer Warm Mission, where the [3.6] — [4.5] > 0.3 
mag color cut combined with the variability criterion will be a solid indication of AGN 
activity. 

The structure function of z > 1 AGNs is well described by a power law with a logarithmic 
slope of 7 = 0.56 ± 0.18 and amplitude S = 0.11 ± 0.02 mag in [3.6] and 7 = 0.44 ± 0.14, 
So = 0.12 ± 0.01 mag in [4.5] for a rest-frame time lag of r = 2 years. These structure 
functions are steeper than those observed in optical bands, and we argue that this is due to 
the inability of either disk or dust emission to respond on the short timescales seen in optical 
studies. The mid-IR structure functions go approximately as r 1 / 2 , which corresp onds to the 



short ti mescale (r <C r mn Hei) be havio r in the damped rando m walk model used by lKelly et al. 



(120091 ). iKozlowski et al.l (120101 ) . and iMacLeod et al.l (|2010| ). If we examine the variability of 
various subsamples, we detect more variability in low luminosity objects observed at shorter 
rest-frame wavelengths. The mid-IR sources corresponding to X-ray and 24 /im AGN have 
structure functions identical to that of z > 1 AGNs. However, this is not the case for the 
counterparts of radio sources which have Sq ~ 0.05 ± 0.01 mag in both [3.6] and [4.5] bands, 
as compared to Sq ~ 0.10 ± 0.01 for either z > 1 objects, X-ray, or MIPS sources. One 
explanation may be that the radio sample has more contamination by non-AGN (e.g., due 
to confusion between radio cores and lobes), another is that the ra dio-selected AGN s are 
accreting in a different mode (lower accretion rates or efficiencies; e.g.. lTasse et al.ll2008l ) and 
that variability in such a mode is inherently lower. 

Expanding the SDWFS survey with new epochs during the Spitzer Warm Mission would 
represent a significant improvement. Given our structure functions, we can estimate the 
results from expanding the SDWFS survey with new epochs. We considered two cases of (1) 
doing one new epoch in 2011 March and (2) doing new epochs in both 2011 and 2012 March. 
The improvements in the mean x 2 are significant, with (x 2 ) = 4 + 3.0S'q/<t 2 and 5 + 4.3S , q/ct 2 , 
respectively, for AGNs, compared to (x 2 ) = 4 or 5 for non-variable objects. These represent 
a significant improvement over the SDWFS baseline, where (x 2 ) = 3 + 1.6Sg/ ' a2 101 AGN and 
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3 for non- variable objects. Figure [TOl shows the increases in completeness for these five and 
six epoch scenarios. Most of the gain for detecting variable sources comes from adding one 
additional epoch. The greatest astrophysical gain from two additional epochs comes from 
being able to better map out and subdivide the structure function by luminosity, redshift, 
and source type to better characterize the physics of AGN variability at these wavelengths. 
In particular, it would confirm the different slope of the mid-IR structure function from that 
observed in the optical. 
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Fig. 1. — Examples of data filtering and image subtraction. All six stamps cover the same 
area (1.2x1 arcmin 2 ) and are centered on one of the most variable AGNs in the Bootes field, 
SDWFS J143032.42+332158.8, with z = 0.43, a 12 = 13.5, and r = 1.0 (see ©. Top row: 
stamp of the original 3.6 fim image (a). Before we do the image subtraction, we filter the 
images in order to convert the Spitzer IRAC's triangle-shaped PSFs into more Gaussian-like 
ones (b). Middle row: subtracted images for the 3.6 /im channel. The AGN brightened 
by about 0.5 mag on a baseline of 4 years between epochs 1 (c) and 3 (d). Black (white) 
color here reflects the source emitting less (more) flux than in the template image defined 
by the average of all four epochs. All non-variable objects are subtracted out. Bottom row: 
subtracted images for the 4.5 /zm channel. The AGN brightened by ~ 0.6 mag between 
epochs 1 (e) and 3 (f), consistent with what is seen in the 3.6 /xm channel. 
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Fig. 2. — Variability in the 3.6 fim (top) and 4.5 /im (bottom) channels as a function of 
magnitude for all objects with four epochs of data shown as contours. The objects are binned 
into 0.1 mag and 0.05 log(f) bins. The contours are drawn for 2, 10, and 100 objects per 
bin, counting from the outer contour. The lower solid line shows the median v m calculated 
in 0.1 mag bins, and the upper three dashed, dash-dotted, and dotted lines show variability 
significances (v — v m )/a = 1, 2, and 3 (from bottom to top) with respect to this median. 
We limit our Photometry Group 2 and 3 samples to objects with [3.6] < 19.7 mag, and 
[4.5] < 19.3 mag (vertical dashed lines). Dots mark ~ 3200 Photometry Group 3 AGNs 
from the modified AGN wedge, selected as either X-ray or 24 /im AGN candidates or having 
a QSO spectroscopic template code. Note that they tend to have higher average variability 
than the typical sources, but with amplitudes easily masked by the photometric uncertainties. 
Open circles mark objects with a u > 10 and r > 0.8. 
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Fig. 3. — Correlation between the 3.6 /im and 4.5 /mi light curves for sources with measured 
[5.8] — [8.0] colors ("Photometry Group 3" sources; see Section [3]) with four epochs of data 
for Variability Level 2, 3, and 4 [an > 2, > 3, and > 4) sources. We only consider sources 
in the shaded, highly correlated r > 0.8 region as true variables, and we use the density of 
the poorly correlated r < 0.5 sources to estimate our false positive rates (Tables [3H5]). 
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Fig. 4. — SDWFS mid-IR color-color diagrams for the Bootes field. The green contours show 
the distribution of all 39,522 Photometry Group 3 sources. The objects are binned in 0.05 
mag bins in both axes, and the contours are drawn at levels of 2, 10, 20, 50, 100, and 150 
objects per bin, counting from the outer contour. Top-left panel: black dashed and solid 
line "wedges" show the original IStern et al.l (120051 ) AGN selection region and the modified 
AGN region we defi ne to encompass the distribution of XBootes X-ray sources discussed in 



Gorjian et al.l ( 120081 ). The black box at color zero is the region occupied by normal stars and 



also nearby galaxies without PAH emission. The blue solid, dotted, and dashed lines show 
the evolution of the mid-IR colors of l ate-type spiral, irreg ular, and early-type galaxies with 
redshift based on the SED models of lAssef et al.l (j2010bl ). The red dotted and solid lines 
show the evolution of a pure AGN and a composite with equal contributions from the AGN 
and its late-type host to the bolometric luminosity. The dots on each track mark redshift 
increments of 0.5 (1.0) from z = to 3 (to 6) for galaxies (AGNs) that start at the z = 0.0 
bulls eye. Other panels: Photometry Group 3 mid-IR counterparts of radio sources (top 
right), X-ray sources (bottom left), and MIPS 24 /xm AGN candidates (bottom right) are 
shown. 
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Fig. 5. — Left panel: mid-IR color-color distribution of all Photometry Group 3, Variability 
Level 2 sources in SDWFS. The conto urs show th e distr ibution of all Photometry Group 3 
sources, as described in Figure HJ The IStern et al.l ( 120051 ) AGN selection region is shown by 
the dashed wedge, and the modified AGN wedge is shown with the solid line. The stellar box 
encompasses sources with (stellar) colors near zero. The majority of mid-IR variable objects 
(76%) are AGNs. Also 11% of objects in the stellar box are variable. Right panel: fraction 
of Photometry Group 2 sources (e.g., sources lacking reliable [5.8] and/or [8.0] photometry) 
that are Variability Level 2 variables as a function of [3.6] — [4.5] color. 
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Fig. 6. — Left panel: mid-IR color-color distribution of the Photometry Group 3, Variability 
Level 2 sources selected as AGNs based on their X-ray (squares), radio (triangles), or 24 
jum emission (dots). The contours and selection regions are the same as in Figure [5j Right 
panel: histograms of the fraction of the Photometry Group 2 X-ray (solid), 24 fim (dashed), 
and radio (dotted) sources identified as Variability Level 2 variables as a function of [3.6] — 
[4.5] color. 
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Fig. 7. — Left panel: mid-IR color-color distribution of the Photometry Group 3, Vari- 
ability Level 2 sources with spectroscopic redshifts from AGES. Stars, circles, and squares 
correspond to objects with stellar, galaxy, or AGN spectroscopic template matches. The 
contours show the distribution of galaxies (gray), stars, and AGNs with AGES redshifts. 
Right panel: fraction of objects with stellar, galaxy, or AGN template matches (dotted, 
dashed, and solid line, respectively) that are the Photometry Group 2, Variability Level 2 
variables as a function of [3.6] — [4.5] color. 
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Fig. 8. — Mid-IR AGN structure functions (Equation (jHJ)) in the observed frame (top) and 
rest-frame (bottom). They are based on ~ 180, 690, and 1100 z > 1 SDWFS objects with 
[3.6] < 16 mag (dots), [3.6] < 17 mag (triangles), and [3.6] < 18 mag (squares), respe ctively. 
We also show the z-band rest-frame structure function from lVanden Berk et al.l ( 120041 ) in the 
lower panels. 
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Fig. 9. — Rest-frame mid-IR AGN structure functions for the objects with [3.6] < 18 mag. 
We divide the objects into two redshift bins, 1 < z < 2 (filled symbols) and 2 < z < 3 
(open symbols), and also into two brightness levels, objects brighter (circles) and fainter 
(squares) than the median for the redshift bin. For time lags of half a year and longer, 
fainter and lower redshift sources show greater variability. The solid lines are the mid-IR 
structure function fits from the bott om panel of Figure El and the dotted lines are the z-band 
rest-frame structure functions from IVanden Berk et al.l (120041 ) . 
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Fig. 10. — Detection completeness for variable AGNs as a function of magnitude for sources 
with high correlation r > 0.8 and either <r 12 > 2 or 4 (Variability Level 2 or 4). We show 
the completeness for the current four epoch SDWFS survey, for the current survey plus one 
epoch in 2011 March (five epochs), and for the current survey plus two epochs in 2011 and 
2012 March (six epochs). 



Table I. SDWFS-DIA Variability Catalog. 



Object Name 


R.A. 


Decl. 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


JVe 


"([3.6]) 


"([4.5]) 


C12 


r 


a 12 


Mask 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 




(mag) 


(mag) 


(mag 2 ) 








SDWFS J143040. 36+322721. 8 


217.668148 


32.456048 


17.666 


17.488 


16.792 


16.070 


4 


0.051 


0.062 


0.001 


0.178 


0.894 


1 


SDWFS J143034. 08+322729.0 


217.642013 


32.458058 


17.937 


17.466 


17.142 


16.928 


4 


0.067 


0.185 


0.002 


0.127 


0.758 


1 


SDWFS J143033. 71+322730. 3 


217.640453 


32.458420 


18.725 


18.352 


17.576 


16.869 


4 


0.181 


0.223 


0.029 


0.718 


0.392 


1 


SDWFS J143047. 11+322729.5 


217.696285 


32.458194 


17.121 


16.995 


16.898 


14.735 


4 


0.036 


0.096 


0.002 


0.511 


0.295 


1 


SDWFS J142951. 47+322726.0 


217.464455 


32.457230 


18.361 


17.875 


17.538 


16.447 


4 


0.135 


0.102 


0.009 


0.627 


0.734 


1 


SDWFS J143015. 91+322733. 8 


217.566284 


32.459402 


19.732 


20.200 


18.262 


99.999 


1 


99.999 


99.999 


99.999 


99.999 


99.999 


1 



Note. — This is an abridged version of the full table, presenting the variability statistics for the combined channels [3.6] and [4.5] only. The electronic table includes additional 
columns showing the variability statistics for the separate channels [3.6] and [4.5], and the photometric uncertainties for the IRAC channels. Here, quantities ii([AT]) and Nc, arc 
respectively, the variability standard deviation in channel [X] and the common number of epochs in both channels used in the calculation of the joint variability. With C"i2i r, 
and (Ti2 we denote the covariance and the correlation between the channels 1 and 2, and the joint variability significance, respectively. The mask is (1) if there is (is no) bright 
star with the 2MASS X-band magnitude within 260(1 — K/ll) arcsec affecting the measurement. The Vega magnitudes presented in the table are measured in 4 arcsec apertures 
(see Ashby et al. 2009). The error code for magnitudes, reflecting no measurement, is 99.999. The error code for v, C\2i r -> an d o\2 ls 99.999. 
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 
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Table 2. The Light Curves for the SDWFS-DIA Variability Catalog. 



Channel 1 Channel 2 



R.A. 


Decl. 


Epoch 1 


Epoch 2 


Epoch 3 


Epoch 4 


Epoch 1 


Epoch 2 


Epoch 3 


Epoch 4 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


217.668148 


32.456048 


17.703 


17.823 


17.741 


17.771 


17.457 


17.518 


17.547 


17.409 


217.642013 


32.458058 


17.894 


18.027 


17.881 


17.952 


17.290 


17.455 


17.566 


17.728 


217.640453 


32.458420 


19.037 


18.793 


18.769 


18.598 


18.721 


18.220 


18.288 


18.371 


217.696285 


32.458194 


17.107 


17.128 


17.132 


17.191 


17.019 


16.995 


17.195 


17.143 


217.464455 


32.457230 


18.385 


18.070 


18.166 


18.148 


17.925 


17.780 


17.696 


17.876 


217.566284 


32.459402 


20.055 


19.751 


20.271 


19.649 


99.999 


99.999 


19.692 


99.999 



Note. — The error code for magnitudes, reflecting no measurement, is 99.999. 
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for 
guidance regarding its form and content.) 



Table 3. Mid-IR Variable Objects in the SDWFS. 



Variability 


Variability Old AGN Wedge Modified AGN Wedge Stellar Box Outside" 


All Abjects a 


Level 


C12 












Number (fraction/false positive 


rate 6 ) of variable objects: Photometry Group 1 




Total 




43468 


76005 45939 144263 


474179 


A 


> 2 


746 (1.7%/28%) 


1015 (1.3%/33%) 555 (1.2%/58%) 832 (0.6%/71%) 


OIL)/ ^X.-L/O/O/ /0) 


Q 
O 


> 3 


397 (0.9%/17%) 


501 (0.7%/22%) 226 (0.5%/54%) 326 (0.2%/56%) 


990R iC\ ^% 1 A7°/r,\ 


A 
I 


> 4 


246 (0.6%/ll%) 


293 (0.4%/17%) 102 (0.2%/51%) 156 (0.1%/53%) 


1071 ( 9% / AA°Z*\ 






Number (fraction/false positive 


rate 6 ) of variable objects: Photometry Group 2 




Total 




29181 


54920 18150 113299 


213594 


2 


> 2 


598 (2.0%/18%) 


774 (1.4%/24%) 133 (0.7%/65%) 512 (Q.5%/74%) 


1557 (Q.7%/50%) 


3 


> 3 


340 (1.2%/9%) 


405 (0.7%/13%) 58 (0.3%/56%) 166 (0.1%/56%) 


668 (0.3%/30%) 


1 


> 4 


219 (0.8%/5%) 


250 (Q.5%/9%) 24 (Q.1%/42%) 71 (0.1%/55%) 


359 (0.2%/24%) 






Number (fraction/false positive 


rate 6 ) of variable objects: Photometry Group 3 




Total 




6677 


12741 7567 19214 


39522 


2 


> 2 


496 (7.4%/6%) 


584 (4.6%/7%) 86 (l.l%/58%) 105 (0.5%/44%) 


775 (2.0%/19%) 


3 


> 3 


310 (4.6%/3%) 


349 (2.7%/4%) 44 (Q.6%/48%) 35 (0.2%/44%) 


428 (1.1%/12%) 


4 


> 4 


204 (3.1%/2%) 


224 (1.8%/3%) 21 (0.3%/33%) 15 (0.1%/30%) 


260 (0.7%/9%) 



Note. — Photometry Groups 1 and 2 are defined in Section [3] only 49% and 86% of objects in these groups have all four IRAC band 
measurements, respectively. The remaining objects lack either [5.8], [8.0], or both measurements, so we are unable to place them on the mid-IR 
color-color diagram (i.e., we were unable to tell if they fall into the modified/old AGN wedge or stellar box). By definition, Photometry Groups 
3 includes only objects with all four IRAC measurements. 

a Column "Outside" includes objects with all four IRAC band measurements only that are both outside of the "Modified AGN Wedge" and the 
"Stellar Box." Column "All Objects" contains all objects, regardless of the number of IRAC band measurements. 

b The fraction of variable sources is a ratio of the number of variable objects above a given level of <ri2 and r > 0.8 to the total number of objects 
in a respective region. The false positive rate is estimated from the number of sources above the same level of <ri2 but with r < 0.5 (see Section|3] 
and Figure [3J. 
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Table 4. Mid-IR Variability of Photometrically-Selected AGNs. 



Variability 
Level 


Variability 
C"12 


MIPS QSOs X-ray FIRST WSRT 


Number (fraction/false positive rate) of variable objects: Photometry Group 1 


Total 
2 
3 
4 


> 2 

> 3 

> 4 


1964 3231 313 2046 
230 (11.7%/5%) 392 (12.1%/4%) 12 (3.8%/8%) 48 (2.3%/13%) 
151 (7.7%/3%) 253 (7.8%/2%) 8 (2.6%/3%) 27 (1.3%/8%) 
107 (5.4%/2%) 164 (5.1%/1%) 4 (1.3%/1%) 17 (0.8%/4%) 

Number (fraction/false positive rate) of variable objects: Photometry Group 2 


Total 
2 
3 
4 


> 2 

> 3 

> 4 


1942 2950 292 1883 
229 (11.896/5%) 391 (13.3%/3%) 12 (4.196/7%) 48 (2.596/10%) 
151 (7.896/396) 253 (8.6%/296) 8 (2.796/3%) 27 (1.4%/696) 
107 (5.5%/2%) 164 (5.6%/l%) 4 (1.4%/1%) 17 (0.9%/3%) 

Number (fraction/false positive rate) of variable objects: Photometry Group 3 


Total 
2 
3 
4 


> 2 

> 3 

> 4 


1637 2217 195 1233 
225 (13.7%/4%) 374 (16.9%/2%) 11 (5.6%/4%) 43 (3.5%/9%) 
151 (9.2%/3%) 245 (11.1%/1%) 8 (4.1%/1%) 25 (2.0%/5%) 
107 (6.5%/l%) 160 (7.2%/l%) 4 (2.1%/1%) 15 (1.2%/4%) 



Note. — The fraction of variable sources is a ratio of the number of variable objects above a given level of a\2 and r > 0.8 to 
the total number of objects in a respective class. The false positive rate is estimated from the number of sources above the same 
level of cti2 but with r < 0.5 (see Section [3] and Figure[3]l- 
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Table 5. Mid-IR Variability by AGES Spectroscopic Classification. 



Variability 
Level 


Variability 

012 


QSOs Stars Galaxies 


Number (fraction/false positive rate) variable objects: Photometry Group 1 


Total 
2 
3 
1 


> 2 

> 3 

> 4 


2950 1121 18161 
334 (11.3%/5%) 20 (1.896/20%) 209 (1.196/25%) 
216 (7.3%/3%) 11 (1.096/7%) 108 (0.6%/15%) 
144 (4.996/2%) 6 (0.5%/5%) 60 (0.3%/10%) 

Number (fraction/false positive rate) of variable objects: Photometry Group 2 


Total 
2 
3 
1 


> 2 

> 3 

> 4 


2820 966 17405 
333 (11.896/4%) 18 (1.996/14%) 200 (1.196/24%) 
216 (7.7%/2%) 10 (1.0%/5%) 102 (0.6%/13%) 
144 (5.1%/2%) 6 (0.6%/4%) 58 (0.3%/8%) 

Number (fraction/false positive rate) of variable objects: Photometry Group 3 


Total 
2 
3 
1 


> 2 

> 3 

> 4 


2192 491 11214 
319 (14.6%/4%) 14 (2.9%/9%) 175 (1.6%/17%) 
211 (9.6%/2%) 10 (2.0%/3%) 93 (0.8%/10%) 
142 (6.5%/l%) 6 (1.2%/3%) 56 (0.5%/5%) 



Note. — The fraction of variable sources is a ratio of the number of variable objects above a given level of u\2 and r > 0.8 
to the total number of objects in a respective class. The false positive rate is estimated from the number of sources above 
the same level of a\i but with r < 0.5 (see Section [3] and Figure l3l. 
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Table 6. Mid-IR AGN Structure Functions 







Structure Function Parameters 


Level 


Channel 


7 


So 


Observed Frame (to = 4 years) 


[3.6] < 16 mag 


[3.6] 


0.42 ±0.14 


0.08 ±0.02 


[3.6] < 17 mag 


[3.6] 


0.52 ±0.16 


0.09 ±0.02 


[3.6] < 18 mag 


[3.6] 


0.52 ±0.14 


0.10 ±0.02 


[3.6] < 16 mag 


[4.5] 


0.28 ±0.11 


0.08 ±0.01 


[3.6] < 17 mag 


[4.5] 


0.40 ±0.10 


0.10 ±0.01 


[3.6] < 18 mag 


[4.5] 


0.45 ±0.14 


0.10 ±0.01 


Rest-Frame (to = 2 years) 


[3.6] < 16 mag 


[3.6] 


0.45 ±0.18 


0.08 ±0.02 


[3.6] < 17 mag 


[3.6] 


0.66 ±0.21 


0.11 ±0.02 


[3.6] < 18 mag 


[3.6] 


0.59 ±0.14 


0.12 ±0.02 


[3.6] < 16 mag 


[4.5] 


0.33 ±0.19 


0.09 ±0.02 


[3.6] < 17 mag 


[4.5] 


0.44 ±0.10 


0.12 ±0.01 


[3.6] < 18 mag 


[4.5] 


0.49 ±0.11 


0.13 ±0.01 


Rest-Frame (to = 2 years) 


1 < z < 2 bright 


[3.6] 


0.45 ±0.17 


0.11 ±0.02 


1 < z < 2 faint 


[3.6] 


0.60 ±0.14 


0.13 ±0.02 


2 < z < 3 bright 


[3.6] 


0.30 ±0.28 


0.09 ±0.02 


2 < z < 3 faint 


[3.6] 


0.19 ±0.11 


0.09 ±0.01 


1 < z < 2 bright 


[4.5] 


0.33 ±0.13 


0.10 ±0.02 


1 < z < 2 faint 


[4.5] 


0.58 ±0.12 


0.15 ±0.01 


2 < z < 3 bright 


[4.5] 


0.13 ±0.15 


0.07 ±0.02 


2 < z < 3 faint 


[4.5] 


0.30 ±0.20 


0.10 ±0.02 
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Table 7. Mid-IR Variability Amplitudes At Fixed 7 = 0.5. 





Amplitude Parameters 


Level 


So [3.6] (mag) 


So [4.5] (mag) 




z > 1 objects (to 


= 4 years) 


14 < [3.6] < 15 


0.07 ±0.01 


0.08 ±0.01 


15 < [3.6] < 16 


0.10 ±0.01 


0.10 ±0.01 


16 < [3.6] < 17 


0.10 ±0.01 


0.10 ±0.01 


17 < [3.6] < 18 


0.10 ±0.02 


0.09 ±0.02 


X-ray 


sources inside the modified AGN wedge (to = 4 years) 


14 < [3.6] < 15 


0.11 ±0.01 


0.10 ±0.01 


15 < [3.6] < 16 


0.12 ±0.01 


0.10 ±0.01 


16 < [3.6] < 17 


0.09 ±0.01 


0.08 ±0.01 


17 < [3.6] < 18 


0.10 ±0.02 


0.11 ±0.03 


MIPS 


sources inside the modified AGN wedge (to = 4 years) 


14 < [3.6] < 15 


0.11 ±0.01 


0.11 ±0.01 


15 < [3.6] < 16 


0.10 ±0.01 


0.10 ±0.01 


16 < [3.6] < 17 


0.09 ±0.01 


0.09 ±0.01 


17 < [3.6] < 18 


0.11 ±0.03 


0.08 ±0.03 


Radio (WSRT+FIRST) sources inside the modified AGN wedge (to = 4 years) 


14 < [3.6] < 15 


0.06 ±0.02 


0.06 ±0.01 


15 < [3.6] < 16 


0.04 ±0.01 


0.06 ±0.02 


16 < [3.6] < 17 


0.05 ±0.01 


0.05 ±0.01 


17 < [3.6] < 18 


0.05 ±0.01 


0.06 ±0.02 



